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ABSTRACT 

We present an analytic approach to predict gas density and temperature profiles in dark matter 
haloes. We assume that the gas density profile traces the dark matter density profile in outer 
parts of the haloes, as suggested by many hydrodynamic simulations. Under this assumption, 
the hydrostatic equilibrium uniquely fixes the two free parameters in the model, the mass- 
temperature normalization and the poly tropic index, that determine the temperature profile. 
This enables us to predict a universal gas profile from any universal dark matter profile. Our 
results show that gas tracing dark matter in the outer parts of the haloes is inconsistent with 
gas being isothermal; on the contrary, it requires temperature to decrease slowly with radius, 
in agreement with observations. We compare our predictions for X-ray surface brightness 
profiles of the haloes and the mass-temperature relation with observations. We find that they 
are generally in a good agreement. We compare the universal profile with the profile and find 
that, although the (3 profile gives a reasonable fit to our predicted profiles, the deviation from it 
can explain many of the observed trends, once we take into account the observational selection 
effects. Our model predicts that the mass-temperature relation does not follow the self-similar 
relation because of the mass-dependent halo concentration. We also predict surface brightness 
profiles of the Sunyaev-Zel'dovich (SZ) effect. We find that fitted to the /3-profile the core 
radii and (3 inferred from the SZ effect are systematically larger than those from the X-ray 
measurement. 
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1 INTRODUCTION 

The structure of the dark matter haloes has been recognized as one 
of the most powerful probes of small scale physics, which may po- 
tentially distinguish between the competing theoretical models. In 
cold dark matter (CDM) models, the structure forms through a se- 
quence of mergers, which gives a characteristic density profile of 
the resulting haloes. 

Numerical simulations led the way in our understanding of 
the structure of the haloes. While there is significant scatter among 
the individual simulated haloes, one finds some similarity among 
them, justifying the introduction of a universal dark matter halo 
profile. The density of such a profile changes steeply in the outer 
parts of haloes, scaling roughly as p oc r~ Q with a ~ 3. The 
slope gradually becomes shallower in the inner parts, converging 
to a value between 1 < a < 3/2 (Navarro, Frenk & White 1996, 
1997; Moore et al. 1998). The transition between the two regimes is 
in general mass dependent, because the lower mass haloes collapse 
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at a higher redshift, when the mean density of the universe was 
higher. This in turn requires a higher central density of the lower 
mass haloes. 

Different dark matter models make different predictions for 
the slope of the inner profile as well as for the value of the central 
density as a function of the halo mass: both CDM and warm dark 
matter (WDM) predict cuspy haloes, but CDM in general predicts 
denser haloes at a given mass and radius. Collisional dark matter, 
however, predicts a central core. 

While the dark matter density profiles are important, they can- 
not be easily observed. In this paper, we focus instead on the gas 
density and temperature profiles, which can be observed through 
the X-ray emission and the Sunyaev-Zel'dovich (SZ) effect. We 
introduce a concept of a universal gas density profile in analogy 
with which an universal dark matter profile describes, in an average 
sense, how the gas density depends on radii. Such a construction is 
as useful as the universal dark matter density profile, in the sense 
that it can parameterize, in a simple way, how different cosmologi- 
cal models differ in their predictions for the gas density profile, and 
place constraints on the models from the observations. 

By far the most common form used for the gas density pro- 
file is the so-called profile, which consists of a central core and 
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a power-law outer slope, (3. While such a profile is useful in obser- 
vational description of the X-ray clusters, it becomes insufficient 
when bolometric SZ observations become available; these observa- 
tions will measure the gas density profile out to the outer parts of 
the clusters, where the predicted slope of the profile changes and 
can no longer be described by a single power law. 

As we discuss in this paper, there are already some indica- 
tions that the j3 profile is insufficient to describe the X-ray data. 
For example, we can explain a correlation between core radii and 
(3, found in observational data, as a consequence of forcing the f3 
profile to describe the underlying gas density profile in which the 
slope increases continuously with radius. 

Another shortcoming of the /3 profile is that it is not derived 
from a dark matter density profile; thus, we cannot relate the ob- 
served X-ray profile to the underlying dark matter profile in simple 
terms. In this paper, we give a prescription of how to compute the 
gas density and temperature profiles from any given dark matter 
density profiles. We show that the prescription is unique (no free 
parameters required) under our assumptions. This property causes 
our model to differ from similar recent constructions by Makino, 
Sasaki & Suto (1998) and Suto, Sasaki & Makino (1998). 

To construct the gas density and temperature profiles, we as- 
sume that the gas is in hydrostatic equilibrium throughout the halo. 
This in itself does not suffice to construct a gas density profile, be- 
cause the dependence of the temperature on radii remains unspec- 
ified. As suming isothermal pr ofile gives a poor fit to the observed 
clusters (Makino et al. 1998) and is not favored observationally 



(Finoguenov, Reiprich & Bohringer 2001). 

Instead, we make an alternative assumption that the gas den- 
sity profile traces the dark matter density profile in the outer parts 
of the halo. While this assumption is qualitatively reasonable, its 
quantitative justification comes from hydrodynamic simulations 
that show remarkable agreement between the two profiles in the 
outer parts of haloes. It is possible that some physics is missing in 
the current generation of numerical simulations, which for the most 
part ignore radiative processes. Gas cooling and star formation can 
affect not only the central regions, but also the global properties of 
the cluster; however, even in that case , the gas profile in the outer 
regions remains relatively unaffected (|Pearce et al. 200d; L ewis et 
al. 2000). 

The remarkable property is that once assuming that the gas 
profile traces the dark matter profile outside the core, and the gas 
obeys a polytropic equation of state, we can predict both the mass- 
temperature normalization and the temperature dependence on radii 
analytically, we find that the predictions agree with observations as 
well as simulations. 

Our assumption that gas temperature is monotonic with den- 
sity limits us to regions outside the core (roughly outside inner 
100 — 200 kpc), because in the inner regions gas temperature is 
often increasing with radius up to 100 — 200 kpc and then mildly 
decreasing in the outer regions. For the most part, however, we 
will compare our results with the observational data for which the 
central part of the cluster has been excised (e.g. Finoguenov et al. 
2001). We are thus justified to make this assumption; consequently, 
our model does not provide constraints on the gas and dark matter 
profiles in the inner parts of the cluster. 

Although the cosmology does not play a major role, for the 
purpose of this paper, we fix the cosmological parameters at the 
best-fittin g values from CMB and large-scale structures (N etterfield 
et al. 2000): fi m = 0.3, Ov = 0.7, Q b = 0.05, h = 0.67, n = 1, 
and as = 0.9. We evaluate all the predictions at the present epoch 
(z = 0). 



The structure of this paper is as follows. In §2, we describe the 
universal dark matter haloes and relations between the virial mass 
and the observed mass. In §3, we introduce a family of self-similar 
gas density and temperature profiles. We show how the assump- 
tions fix free parameters in the model. After discussing relations 
between the virial temperature and the observed gas temperature, 
we predict gas density and temperature profiles of haloes explic- 
itly. In §4, we compare the X-ray surface brightness profiles de- 
rived from the predicted gas profiles with the (3 profile. We study 
the observational selection effect on the comparison in detail. In 
§5, we discuss the surface brightness profiles of the SZ effect. In 
§6, we predict the mass-temperature relation, and compare it with 
observations. Finally, we summarize our results in §7. 



2 UNIVERSAL DARK MATTER DENSITY PROFILES 
2.1 Basic properties 

Many high-resolution TV-body simulations suggest that the dark 
matter density profile, pdm(r), is well described by a self-similar 
form. The profile is expressed with a characteristic radius r s : 



Pdm{r) = p B ydm(r/r B ), 



(1) 



where p s is a normalization factor in units of mass density, while 
ydm(x) is a non-dimensional function representing the profile. We 
have chosen p s so as to represent a characteristic density at a char- 
acteristic radius, r — r s . The characteristic radius describes a typi- 
cal scale at which the profile slope changes from the outer value to 
the inner value. Often it is close to the radius at which the slope is 
a = 2 (Navarro et al. 1996, 1997; Moore et al. 1998). 

Since the dark matter density profile is self-similar, the dark 
matter mass profile is also self-similar. The dark matter mass en- 
closed within a radius r is 



M(< r) = Aivp s r s m (r/r s ) , 

where m(x) is a non-dimensional mass profile given by 



m(x) 



duu 3/dm(w) 



(2) 



(3) 



(/dm(ii) must be shallower than u -3 in the central region in order 
for the mass to converge. 

We define total dark matter mass, M v i r , as the mass within the 
virial radius, r v i r , i.e., M v i r = M(< c), where 



' vir 

r B 



(4) 



is a dimensionless parameter called the concentration parameter. 
Evaluating equation (Q) at the virial radius, we fix the normalization 
factor, p s , as 



ps = C 



47rr 3 ir m(c) ' 



(5) 



We calculate t he virial radiu s, r v i r (Af v ir, z), with the spheri- 
cal collapse model ( Peebles 198C ), 

1/3 

(6) 



M vir 



(47r/3)A c ( 2 )p c (z) 

where A c (z) is a spherical overdensity of the virialized halo within 
r v ir at 2, in units of the critical density of the universe at z, p c {z). 
In the Einstein-de Sitter universe model, A c (z) is independent of 
z, and is exactly 187r 2 ~ 178. Lacey & Cole (1993) and Nakamura 
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& Suto (1997) provide A c (z) for arbitrary cosmological models. 
In our cosmological model with Sl m = 0.3 and Q,a = 0.7, the 
spherical overdensity of the virialized halo at the present is 100: 
A c (z = 0) = 100. 

The only parameter in our framework is the concentration 
parameter, c. Many iV-body simulations show that c decreases 
gradually with the virial mass. Assuming the mass dependence of 
c known, for example, from iV-body simulations (Bullock et al. 
2001 ; Eke, Navarro & St einmetz 2000 ), or from the non-linear dark 
matter power spectrum (Seljak 200C), one can predict pdm(r) f° r 
any Mvir, any z, and in principle any cosmological models. This 
is the main advantage of the universal dark matter density profile 
paradigm, and explains its popularity in recent work. Although the 
universal profile does not explain the dark matter profiles o n the 
object-to-object basis ( Ting & Suto 2000} Klypin et al. 2001), it is 
still successful in fitting the averaged form of the dark matter den- 
sity profiles with reasonable accuracy. 

So far, we have discussed the functional dependence of the 
profile on the characteristic radius. For a complete description, 
however, we need the functional form of ydm(x). We adopt 



X a (l + x) 



(7) 



The asymptotic profile in the x S> 1 regime is ydm{x ^> 1) = 
x~ 3 . This is the most common form found in the iV-body simula- 
tions, altfimigJ^C£rtter_arc™^ et 
al. 2001). In the i < 1 regime, we have ydm(x <C 1) = x~ a . Par- 
ticularly, a = 1 corresponds to the profile proposed by Navarro et 
al. (1996, 1997), while a = 3/2 corresp onds t o the one proposed 
by Moore et al. (jl998b and Jing & Suto (kOOOh. We will focus on 



these two particular cases throughout this pape r. 

Using these parameters, Suto et al. (1998) have evaluated the 
relevant integrals for a 



1: 



m(x) 



du 



■- ln(l + x) 
m(u) _ 



1 + x' 
ln(l + x) 



and for a = 3/2: 



m(x) = 2 In (VS + s/TTx) - 2 \J^~ 



(8) 
(9) 

(10) 



We use J Q duu~ 2 m(u) in next section. 

There are several different empirical fitting formulae for the 
concentrat ion pa rameter in the literature. A recent compil ation is 
Eke et al. (EoOOb. For a = 1, we will use the one of Seljak (feOOOh, 



Mvir 



© 



-1/5 



(12) 



For q _ 3/2, we reduce the concent ration parameter by a factor of 
1.7 (|Seljak 200$ |Klypin et al. 200l|). 



2.2 Overdensity radius and mass 

X-ray observations are not sufficiently sensitive to measure X-ray 
surface brightness and gas temperature profiles out to the virial ra- 
dius. Therefore, instead of measuring the virial mass, many authors 



measure the mass within an overdensity radius, rg, at which the 
dark matter density is 5 times the critical density of the universe: 



rs 



M(< r s ) 



(4tv/3)S Pc (z) 



1/3 



(13) 



Note that the virial radius at z — corresponds to rioo in our 
cosmological model, as A c (0) = 100 in equation (^). Since we 
assume the dark matter density profile, and hence the mass profile, 
known, we can relate M v i r to M(< rs) by solving the equation 



M(< r s ) = 



m(rs/r s ) 


M vir = 


m(c) 





m(crs/r viv ) 
m(c) 



M vir 



(14) 



Figure [l] plots M(< rsoo) as a function of M vlr for a = 1 
and 3/2. The mass dependence arises solely from the mass de- 
pendence of the concentration parameter given by equation (p^|). 
Since we will primarily be looking at haloes with the mass be- 

Mq, the relation between the 



tween 10 A h 



" 1 M and 10 15 /i" 



two is roughly M(< rsoo) ~ M v i T /2. The mass measurement 
with the gas density profiles out to r 5 oo is thought to be accurate 
(e.g. Evrard, Metzler & Navarro 1996); this is why this radius has 
been frequently used in the literature. Therefore, we will use rsoo 
and M(< rsoo) for comparison of our theoretical predictions with 
observations. 

Figure [l] also plots a ratio of r$oo to the virial radius as a func- 
tion of A/vir. One finds that rsoo is half the virial radius. Note that 
there is some ambiguity in the definition of the virial radius. Some 
authors use a fixed overdensity, 8, regardless of cosmology, and 
others use the spherical collapse value, A c , that we use. 

Although we will continue to use the concentration parame- 
ter defined by equation (Q), we may also introduce a more model- 
independent quantity like C500 = rsoo/r s , which does not depend 
on the definition of the virial radius. For the cosmology used here, 
C500 is half the concentration parameter. 



3 SELF-SIMILAR GAS DENSITY AND TEMPERATURE 
PROFILES 

3.1 Polytropic gas model in hydrostatic equilibrium 

Since we assume the dark matter density profile self-similar, a 
gas density profile, p gas (r), would also be self-similar, if no ad- 
ditional scale is introduced (this translates into a requirement that 
the density-temperature relation is a power law). Hence, we have 



(r) 



(0)2/ g as (r/r s ), 



(15) 



where p gas (0) is the gas density at r = 0. We parameterize the 
gas density profile in this way because, in contrast to the dark mat- 
ter density profile, a flat gas density core near the centre of haloes 
is well observed, and is expected if t he ga s temperature is not a 
strong function of radii. Makino et al. ( 1998) and Suto et al. ( 1998) 



show that the hydrostatic equilibrium between the gas pressure and 
the self-similar dark matter potential gives the self-similar density 
profile. 

To take phenomenologically into account the effect of the gas 
temperature gradient which is found both in observations and in 
hydrodynamic simulations, Suto et al. ( 1998 ) parameterize the gas 



pressure in a polytropic form, P gas oc p gas T g: 
rameterization gives a self-similar temperature profile, 

T gas (r/r s ) =T gas (0)^" 1 (r/r s ). 

Then, the hydrostatic equilibrium equation, 



oc p^as- This pa- 



(16) 
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Figure 1. Ratio of the virial mass, M v ; r , to the overdensity mass at 7*500, 
M500 = M(< 7*500), as a function of M v i r . Ratio of 7*500 to the virial 
radius, r v - lr , is also shown. The solid lines represent a = 1, while the 
dashed lines represent a = 3/2. 



-1 ^gi 



c/r 



= -G 



M(< r) 



gives the differential equation for the gas profile, y sas (x), 



dy^s (r/r s ) 
dr 



7 — 1 \ Gfim p M viT 



fc B T gas (0)r 2 



i(r/r 



m(c) 



(17) 



(18) 



Here, M v \r is the virial mass, and m{x) is the dimensionless 
mass profile given byequation (^). Since 7/ gas (0) = 1 (see equa- 
tion ( |l5T)), equation (|l8j) is formally solved for y eas {x) as (S uto et 
al. 1998) 



7 
Vg&s 



7-1 



m(c) 



du 



777(77) 



where 



Gjj,m p M v i r 
3r v irfeBT 1 gas (a;) 



(19) 



(20) 



is the normalization factor relating T gas (x) to M v i r . Note that 

The hydrostatic equation ( |l 8[ ) does not include p gas (0), but 
77(0) representing T gas (0); thus, 77(0) is one of the two free param- 
eters in the model. The other parameter is the polytropic index, 7. 
As we will show below, by requiring the gas and the dark matter 
profiles to agree outside the core, we can determine both the pa- 
rameters. 



3.2 Gas tracing dark matter outside the core 

We assume that the gas density profile traces the dark matter den- 
sity profile in the outer region of the halo. This assumption is the 
main ingredient of our model presented in this paper. 

Qualitatively, this assumption is reasonable, as the gas and the 
dark matter obey similar equations of motion to each other. The 



main difference is that the gas falling into the cluster potential well 
is shocked and then heated, while the dark matter, being collision- 
less, does not undergo shock heating. What it does undergo is orbit 
mixing, which causes multiple streams to appear at the same physi- 
cal position. This gives the velocity dispersion, an effective temper- 
ature of the dark matter. Taking moments of the collisionless Boltz- 
man n equation leads to the Jeans equations (B inney & Tremaine 
1987) which are formally similar to the gas fluid equations. The 
main difference is that, for the collisionless system, the effective 
pressure is the velocity dispersion tensor that can be anisotropic, 
while the gas pressure is isotropic. 

An even better reason for our assumption is that many hy- 
drodynamic simulations have observed the gas density profile trac- 
ing the dark matter density profile outside the core, r <; r v i r /2 
(Navarro, Frenk & White 1995; Bryan & Norman 1998; Eke, 
Navarro & Frenk 1998; Frenk et al. 1999; Pearce et al. 2000; Lewis 
et al. 2000; Yoshikawa, Jing & Suto 2000 ). For example, the Santa 
Barbara cluster project (Frenk et al. 1995 ), the comparison between 
a dozen numerical codes for a single cluster, shows this property 
with 10% or better accuracy. Other larger samples also show the 
agreement. While additional physics missing in the non-radiative 
simulations, such as cooling and star formation, may change the 
structure of the gas, the outer regions appear to remain relatively 
unaffected (Pearce et al. 200C ; Lewis et al. 2000). 

We require the gas density profile, y gas (s), to match the dark 
matter density profile, ydm(x), at a certain matching point, a**. We 
impose this by requiring slopes of these two profiles to be the same 
at x*. 



dlny d m{x) 



In a; 



dlny gss (x) 



In a; 



(21) 



where s* denotes an effective slope of the dark matter density pro- 
file at x#. Substituting ydm(x) for equation (^), we obtain the ef- 
fective slope at 1,, 



a + (3 - a) 



1 + x, 



(22) 



This constraint, equation (pl|), fixes one of the two parameters of 
the model, 7*(0), the normalization factor of the mass-temperature 
relation given by equation (po[). 

Substituting y gSlS (x) in equation for equation (|l9|), we 
obtain a solution for the mass-temperature normalization factor at 
the centre, rj(0), as 



77(0) 



■{(?) 

+3(7- 



~ 1 m(x„) 
~ 1 m(c) 



In the isothermal limit, 7 



m(c) 
1, we have 



dii 



777(71) 



r?(0) 



x, 1 m(x f ) 



7( C ) 



(23) 



(24) 



Note that this gives 77(0) = 1 for x* = c in the limit of large 
concentration, c — *• 00. 

In order for the two profiles to agree over a wide range, the 
two slopes should agree not only at x = x t , but everywhere in the 
outer region of haloes, so within a factor of two above and below 
the virial radius (c/2 < x < 2c) according to the simulations. 
In this way, the solution for 77(0) will not depend on x*. This re- 
quirement will not be satisfied in general; however, we still have 
one free parameter left, the polytropic index, 7. We fix it by requir- 
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Figure 2. Mass-temperature normalization factor at the centre, r)(0), pre- 
dicted by the polytropic gas model with a polytropic index of 7. r](0) is 
plotted as a function of the matching point, x*, in units of the concen- 
tration parameter, c. The solid lines represent a = 1, while the dashed 
lines represent a = 3/2. From left to right panels, M v i r = 10 15 , 10 14 , 
and 10 13 Mq are plotted. Each mass corresponds to the concentration pa- 
rameter of cnfw = 4.1, 6.5, and 10.4 for a = 1, respectively, while 
c = c NFW /1.7for a = 3/2. In each panel, 7 = 1.0, 1.05, 1.1, 1.15, 1.2 
and 1.3 are plotted from bottom to top lines. 



ing the solution for 77(0) not to depend oni„ Our model then no 
longer contains any free parameters. Using the solution for r/(0), 
we calculate the gas temperature at an arbitrary radial point with 
r)(r)=r } (0)y^ 1 (r). 

Figure plots 77(0) as a function of the matching point, a;*, 
for a variety of polytropic indices, 7. It follows from this figure 
that 7 = 1.1 — 1.2 satisfies the requirement that 77(0) is indepen- 
dent of x*. The isothermal case, 7 = 1, fails this requirement; the 
isothermal gas density profile cannot be made similar to the dark 
matter density profile in the outer parts of the halo. 

The best-fitting value of 7 depends on c very weakly, changing 
from 1.1 to 1.2 as c varies from 4 to 10 (or the mass from M v j r = 
1O 15 M to W 13 M Q ). We find a simple linear fit for 7 = 7(c), 

7 = 1.15 + 0.01 (c N fw - 6.5) , (25) 

where cnfw denotes the concentration parameter for the dark mat- 
ter density profile with a = 1. Equation ( p5| ) gives an adequate 
value to both a = 1 and a — 3/2 profiles. Hence, 7 is no longer 
free, but fixed. Here, note that the concentration for a — 3/2 is 
given by c = cnfw/1-7. 

Furthermore, combining equation (^) with (0) gives 77(0) as 
a function of the concentration parameter. Figure pfplots 77(0) as a 
function of cnfw- We find that the dependence of 77(0) on cnfw 
is well fitted by 

77(0) = 0.00676 (cnfw - 6.5) 2 +0.206 (cnfw - 6.5)+2.48,(26) 



for a = 1, and 



Figure 3. Mass-temperature normalization factor at the centre, r/(0), as a 
function of the concentration parameter for a = 1, cnfw- Th e solid line 
represents a = 1, while the dashed line represents a = 3/2. The dotted 
lines are the best-fitting lines given by equations ( p6| and (jz^). 



for q = 3/2; thus, the higher concentration implies the central 
gas temperature being higher than the gas temperature at the virial 
radius. Since by construction the particular choice of x, does not 
affect our analysis, we take x* — c (or r* = r v i r ) for definiteness 
in the following. 

Figure ^ shows histograms of observationally measured 
polytropic indices . The top panel plots the whole sample from 
Finoguenov et al. (2001). There is a peak in the population around 
7 ~ 1.1. The mean value is 7 = 1.15. This mean value is close 
to our best-fitting value given by equation (|25|). The bottom panel 
of this figure plots the hotter objects than 3 keV in the solid line, 
and the cooler ones below it in the dashed line. The distribution is 
basically the same for these different temperatures; thus, there is no 
significant dependence of the polytropic index on the gas tempera- 
ture. This also agrees with our model. 



3.3 Mean gas temperature of haloes 

Since X-ray observations do not measure the gas temperature itself 
but the emission-weighted mean temperature, we use the latter as a 
representative gas temperature of the halo. If the halo is isothermal, 
then the gas temperature and the emission-weighted mean temper- 
ature are the same; however, the temperature gradient makes them 
different. We predict the emission- weighted mean temperatur e, Tx, 
with this spherical polytropic gas model as (Suto et al. 1998) 



77(0) = 0.00776 (cnfw - 6.5) 2 +0.264 (cnfw - 6.5)+3.07,(27) 



J d^pg as A(T gas )r gas 
/ dV Pgas A(T gas ) 



T g as(0) 



J Xmax x 2 dx[y saB (x)} 



(7+3)/2 



(28) 
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Figure 4. Meas ured distribution of the polytropic index 
Finoguenov et al. (2001). The top panel plots the whole 
different bins: the solid line is A7 = 0.1 bin, and the dashed 
A7 = 0.05 bin. The mean value of 7 is 1.15. The bottom 



sam jli 



the temperature dependence of 7. The solid line shows the 
3 keV, while the dashed line shows the cooler ones below 3 
no clear evidence that the polytropic index changes with temper iture 



7, from 
in two 
line is 
lanel plots 
clusters above 
There is 



i keS. 



where we have approximated the cooling function, A(T 
the bolometric bremsstrahlung, A(T , gas ) oc Tgls oc y g Zs 
that p gas oc y gas , and T gas oc j/^al 1 • Since we d( 
diet the X-ray luminosity but the X-ray temperature, thi; 
mation for the cooling function seems sufficient. 

For a practical evaluation of Tx, we should care how 
the integration boundary, a; max . As y gEls (x) oc x s * ~ x~ 



Sas; , with 
)/2 .Note 
attenjtpt to pre- 
approxi- 



l(7+3)/2 



-3( 7 +l)/2 



to choose 
for x 2> 

a; -3 ; thus, the integral in 



1, x [y gas (x)Y' oc a; 

the denominator of equation ( |28| ) converges rapidly as 
numerator converges even more rapidly than the 
long as 7 > 1. This rapid convergence implies that a 
of x ma x has no systematic effect on Tx. Therefore, we ca 
compare Tx calculated from our model with observations 
simulations. The bottom panel of figure 5 plots the sensiti' 
to :r max . The plotted quantity, rjx, is the emission-weigi 
temperature normalization factor defined by equation ( 
and is equivalent to Tx. We find that Tx is quite insensith e 
for aimax k, 0.2c ~ 0.4r 50 o/r s . We use ir max = c for 

To compare our predictions with observations an 
tions, we replace the normalization factor of the mass- 
relation, 7/(r) (equation (po|)), with the emission-weigl 
normalization factor, 



J'/x 



J dy P g as A(T gas ) ?? (r) 
/ d^pI as A(T gas ) 

~f*—xZdx[y gas (x)f +1)/2 



7,(0) 



f*" 1 ™ x 2 dx[y gas (x)] 



(7+3)/2 



Several hyd rodynamic simulat i ons have calculated this 
tion factor (Evrard et al. 1996; Bryan & Norman 1998 



The 

as 

c choice 
robustly 
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mass- 
below, 

tO X ma x 

de: initeness. 

1 simula- 
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denon linator 



t speci fl 



h ed 



(29) 

n Drmaliza- 
Eke et al. 



Figure 5. 

of the 
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panel is 
profile for 
top to 



Predicted mass-temperature normalization factors as a function 
upptr boundary of the integral, 

Srnaxi in equations <§§) and (0). The 
pa lei is the emission-weighted normalization factor, while the top 
thp mass-weighted one. The solid lines use the dark matter density 
t = 1, while the dashed lines use the one for a = 3/2. From 
hottbm lines, M = 10 15 , 10 14 , and 10 13 Mq are shown. 
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(3/21/t 
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Sine 

cantly weighted 
to the nu: nerical 
significai it 
tions, in 
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dispersion of 7/x among different hydrodynamic simula- 
ontra st to the mass-weighted mean normalizatio n factor, 
defirjed as ( Yoshikawa et al. 2000| ; Thomas et al. 2001) 

J*m*x x 2 dx yl^x) - 



J dVp gsl sV(r) 



This is 
emission 
observed 
the nume - ical 
weighted 
converge 
7 > 1, 
ity of r\_ 
we use 



, but 



(2001) h; .ve 



hikawa et al. 200C ; Thomas et al. 2001 ). For exam ple, rtx 
s to the emission-weighted 7 in Yoshikawa et al. ( |200o| , 



in Bryan & Norman (1993), and (3/2)/3 TM in Eke et al. 



e the emission-weighted mean temperature is signifi- 
toward the central region of haloes, it is sensitive 
resolution ( Lewis et al. 2000 ). We find that there is 



/ dV Pg 



77(0) 



J*™« x 2 dx y gas ( x ) 



(30) 



It ss 



weighted towards the central region of haloes than the 
weighted normalization. Unfortunately, 7y m is not directly 
from the X-ray observations, while it is less sensitive to 
resolution of simulations. In contrast to the emission- 
normalization, the denominator of equation ( |30[ ) does not 
but diverges logarithmically. The numerator converges if 
very slowly. The top panel of figure 5 plots the sensitiv- 
to x max . We find that 77 m is rather sensitive to a; max ; thus, 
tlie virial radius for the integration boundary: x max = c, 
tocomp; re the predicted rjm with simulations, for Thomas et al. 
calculated rj m within the virial radius. 
Usii ig equation ( p9[ ) and (^C ), we calculate 77X and 77m as a 
function jf the halo mass. Figure 61 compares the predictions with 
the simu] ations for the dark matter profiles with a — 1 and 3/2. To 
facilitate the comparison, we plot the dependence against the virial 
mass rather than the concentration parameter, but the two are re- 
lated to t ach other through equation (|l2|). We find that 7fx > 77 m , 
i.e., Tx p T m , as observed in hydrodynamic simulations carried 
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out by Yoshikawa et al. (200C) and Thomas et al. (2001). This is 



because the emission-weighted mean temperature is much more 
weighted towards the central region than the mass-weighted one. 
The gas temperature in the central region is higher than in the outer 
region, as long as 7 > 1. 

We also find that the normalization factors increase as the 
mass decreases, more so for rjx than for 77m . This is caused by 
the concentration dependence of the halo mass. The gas density in 
the core for higher c (smaller mass) is larger than that for lower c 
(larger mass). This in turn gives a smaller mass halo a higher tem- 
perature because of the assumed polytropic relation between the 
gas density and temperature. 

Many hydrodynamic simulations in the literature give the 
mass-temperature normalization factors weighted by either the 
emission or the mass. Figure ^ compares our predictions with those 
obtained from the simulations. Overall, the predicted normaliza- 
tions agree with the simulations in t he hi gh mass regime (except 
for simulations by Bryan & Norman (199£ ) who predict 20% lower 
normalizations compared with the other simulations and with our 
results). 

Our predicted emission-weighted normalizations are higher 
than the simulated ones for smaller mass haloes. It is possible that 
a poorer numerical resolution for smaller mass haloes causes the 
disagreement. It is also possible that the mass dependence has not 
been noticed, as most of previous work has relied on a self-similar 
model for comparison, in which there is no mass dependence of the 
concentration parameter. The predicted mass-weighted normaliza- 
tions which are less sensitive to the resolution issues agree better 
with the simulations. We present the comparison with observations 
in the next section. 



3.4 Predicted gas density and gas temperature profiles 

We plot the dark matter (dashed lines) and the gas density (solid 
lines) profiles in figure [^for M = 10 15 , 10 14 , and 1O 13 M from 
left to right panels, a = 3/2 is plotted in the top panels, while 
a = 1 is plotted in the bottom panels. By construction, the gas and 
the dark matter density profiles are very similar to each other in the 
outer parts of the haloes. 

The gas density profile is always shallower than the dark mat- 
ter density profile in the central part of the halo, developing an ap- 
proximate core. This is because the large thermal gas pressure in 
the centre balances the force of gravity exerted on the gas fluid 
element even for a nearly uniform gas density. The predicted gas 
profiles agree with numerical simulations v ery well: the sim ulated 
gas profil es are also smoother in the centre dEke et al. 1998|; Frenk 
et al. 1999TFearce et al. 2000t lewis et al. 2000|; Y oshikawa et al. 
2000). 

Figure H plots the gas temperature profiles for M = 10 15 , 
10 14 , and 10 3 Mq from top to bottom lines. The left panel shows 
q = 1, while the right panel shows a — 3/2. The predicted tem- 
perature pr ofiles also agree with the simulated te mperature profiles 
very well ([Mavarro et al. 1995|; |Eke et al. 19981; Bryan & Norman 



1998; |Frenk et al. 1999| ; Pearce et al. 2000} |Yoshikawa et al. 2000[ 
Thomas et al. 200 1[ ). 

The predicted gas density profile for M — 10 15 Mq isre- 
markably similar to the simulated cluster by Frenk et al. (1999) 
which has M = 1.1 x 10 Mq. The predicted temperature profile 
also agrees with the simulated one; both show the gas temperature 
declining with radius, which is described b y 7 — 1.1 — 1.2. Usinga 



1.6 



n 1 111 




10 14 I ( 

M vir [h-»M a ] 



Figure 6. Predicted mass-temperature normalization factors as a function 
of the virial mass, M v - lr . The emission-weighted normalization factor, rfx 
(equation (j^)), is plotted in the bottom panel, while the mass-weighted 
one, r; m (equation (po|)), is plotted in the top panel. The thick solid lines 
represent a = 1, while the thick dashed lines represent a = 3/2. The thin 
solid lines indicate the normalization fa ctors d erived from hydrodynamic 
simulations carried out by Thomas et al. (2001 1 (upper li ne in t he top panel 



and lower line in the bottom panel) and Yoshikawa et al. ( 2000 ) (lower line 
in the top panel and upper line in the bottom panel). 



The simulations discussed above do not include radiative pro- 
cesses. One may worry that including cooling and star formation 
may c hange the results and make the agreement worse. Lewis et al. 
(2000) have investigated this point, and conclude that the addition 



larger sample simulated by Thomas et al. (2001 1, we find the best- 
fitting value around 7 = 1.1, which agrees with our prediction. 



of these processes modifies not only the central properties, but also 
the global structure of the haloes to some extent. Nevertheless, they 
find that 7 = 1.1—1.2 gives good fit to their simulated results out- 
side the core and inside the virial radius, regardless of the inclusion 
of the cooling effect. 

We conclude that our predicted profiles with the polytropic 
index fixed by equation (E5j) fit both the observations and hydro- 
dynamic simulations very well, outside the core. Recent Chandra 
(McNamara et al. 2000; Allen, Ettori & Fabian 2001; Allen et al. 
2001) and XMM-Newton (Arnaud et al. 2001a, 2001b; Kaastra et 
al. 2001) observations indicate that inside the core the gas tempera- 
ture rises with radii. Our model with 7 ~ 1.15 cannot describe this; 
thus, we will not attempt to extend our model to the region inside 
the core. 



4 COMPARISON WITH OBSERVED X-RAY SURFACE 
BRIGHTNESS PROFILES 

4.1 Spherical /3 profile 

Historically, X-ray observations have been interpreted with the 
spherical /3 profile, 
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Figure 7. Predicted gas density profiles (solid lines) shown with univer- 
sal dark matter density profiles (dashed lines). The top panels use the dark 
matter density profile for a = 3/2, while the bottom panels use the one for 
a = 1. From left to right panels, each panel shows the case of M = 10 15 , 
10 14 , and 10 13 Mq, respectively. The corresponding concentration param- 
eter, c, to the mass is quoted. 



Figure 9. Predicted X-ray surface brightness profiles (solid lines), and the 
best-fitting spherical fj profiles (dashed lines). The fit is performed out to 
7*500- Note that rsoo is half the virial radius (see figure Mj). The top panels 
use the dark matter density profile for a = 3/2, while the bottom panels 
use the one for a = 1. From left to right panels, each panel shows the case 
of M = 10 15 , 10 14 , and 10 13 M Q , respectively. 
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Figure 8. Predicted gas temperature profiles. The left panel uses the dark 
matter density profile for a = 1, while the right panel uses the one for 
a = 3/2. From top to bottom lines, M = 10 15 , 10 14 , and 1O 13 M 
are shown. The corresponding concentration parameter, c, to the mass is 
quoted. 



/x(r) = Jx(0) 



1/2-3% 



(31) 



where r c is the projected X-ray core radius. This profile had a re- 
markable success in fitting observed X-ray surface brightness pro- 
files, except for the very central parts of haloes, where the radiative 



cooling effect seems to play a major role. Observers often exclude 
this region, when they fit the measured X-ray profile to the (3 pro- 
file. 

To compare our predicted profiles with observationally mea- 
sured X-ray profiles, we fit the core radius and /3x to our predicted 
X-ray surface brightness profile, 



Ix[x) oc 



dl Pg as A(T gas ) 

1 

dl t/ g 



+ P 



(7+3)/2 



(32) 



Figure ^ plots the predicted X-ray profiles, and the best-fitting j3 
profiles. We find that the predicted X-ray profile i s quite similar to 
the P profile ( |Makino et al. 199S| ; |Suto et al. 199^ ). 



4.2 Examining observational selection effects 

The main issue in comparing quantitatively the predicted X-ray 
profiles with the observed ones is out to which radius the observed 
X-ray profile is typically fitted to the (3 profile. Our predicted pro- 
file has a slope decreasing from to —3 continuously with radius, 
corresponding to /3x varying from to 1 ; thus, the larger outermost 
radius we use for the fit, the larger /3x we obtain. 

The top panel of figure [l(j shows how sensitive fix is to the 
maximum radius, r max , that is used for th e fit. We find that fix 
increases with r max . Navarro et al. (1995) also observe this ef- 



fect in their hydrodynamic simulations. We find that the sensitiv- 
ity of fix to r max is comparable to the mass dependence of fix- In 
other words, we can create any correlation between fix and mass 
by changing r max slightly but systematically. The bottom panel of 
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Figure 10. Sensitivity of the parameters of the spherical /3 profile to the 
maximum radius, r max , used for fitting the X-ray surface brightness profile 
to the /3 profile. The maximum radius is in units of the virial radius. The 
top panel shows the X-ray outer slope, /3x> while the bottom panel shows 
the X-ray core radius, r£. The solid lines, dashed lines, and dotted lines 
indicate the virial mass of f0 15 , f0 14 , and 10 13 Mq, respectively. The 
corresponding concentration parameter, c, to the mass is quoted. The dark 
matter density profile for a = 1 is used. 



Figure 11. Observational maximum radii, r max , used for fitting the X-ray 
surface brightness profiles to the (3 profile, as a function of the emission- 
weighted mean temperature, Tx . rmax is in units of rsoo ■ The circles are 



the data taken fr om M ohr et al. (19991, while the triangles are from Hels- 
don & Ponman ( [?00o| ). The solid line plots an estimate of the temperature 
dependence of the selection effect on r max caused by the background noise 
(equation (^)). Note that rsoo is half the virial radius (see figure [jj). 



the figure plots the sensitivity of the core radius to r max . We find 
that r* also increases with r max . 

To make a meaningful comparison between our predicted X- 
ray profiles and observations, we have to characterize the observa- 
tional selection effect on r max . In X-ray observations, r max is the 
radius at which the background noise starts to dominate the X-ray 
surface brightness. Hence, we roughly estimate r max by equating 
■fx(fmax) to the background noise level, BG, 



BG = I: 



x(r D 



\ 2 /2 / 'max \ 

« p gas (0)T gas — J 



(33) 



where s* ~ —3 is a slope of the gas density profile given by equa- 
tion ( p2[ ) with x, = c. As c oc Tg^ 10 , p gas (0) oc c 3 oc Tg^ 10 , 
and r s = r v i r /c oc Tg/s 5 - After some algebra, we obtain r max in 
units of the virial radius, r v i r , 

r max rT n3/10+l/f4s»+2) / n^y\l/(2s*+l) 



ocT 3/10+l/(4 S .+2) (jBG) l/ 



(34) 



If we assume BG not to vary significantly from observation to ob- 
servation, we obtain r max /r v i 



n 3/10+l/(4s„+2) 
1 gas 



T, 



1/5 



Figure plots r max in units of r 50 o for observed clusters 
and gr oups. We take the data from Mohr, Mathiesen & Evrard 
( |1999| ) and Helsdon & Ponman < {200C| ). We calculate r 50 o from the 
emission- w eighted mean temperature, using the normalization of 
Mohr et al. jl999j ), 

r50 ° = 2 ' 37 ( i^v^ V2 1 Mpc - (35) 



Despite the large scatter, there is some evidence that 
k/?"5oo increases with Tx. The figure also plots (solid line) an 



estimated scaling of the selection effect on r max with Tx (Eq.(|34J)). 
We have normalized it so as to have r max = 0.8rsoo at 6 keV, 
which agrees with the observational data. Equation ( p4[ ) seems to 
capture a rough tendency of the actual selection effect, while the 
large scatter around the solid line in the figure would be due to the 
background noise level, BG, varying substantially from observa- 
tion to observation (or clusters being not universal on the individual 
basis). 

We summarize the selection effect as follows: For the cooler 
haloes, the smaller r max /r5oo is used for fitting the observed X- 
ray surface brightness profiles to the (3 profile. This effect leads 
us to have smaller /3x and smaller core radii for the cooler haloes 
than we have for the hotter haloes, according to our predictions. We 
see below more quantitatively how the selection effect affects the 
comparison. 



4.3 X-ray core radii 

Figure |^| compares our predicted X -ray c ore radii, r*, with the 
observational data from Mohr et al. ( 1999[ ) and Finoguenov et al. 
(2001). We have chosen the samples from Mohr et al. ( 1999[ ) to 
which the double-/? analysis is not applied. The thick lines plot 
the predictions corrected for the selection effect due to the back- 
ground noise (equation (Q)), while the thin lines plot the predic- 
tions not corrected. For the latter, we have used r max = rsoo re- 
gardless of the temperature. The selection-effect correction makes 
W smaller, and thus makes the core radii smaller than those with- 
out the correction, as the gas density profile is shallower in the inner 
region than in the outer region (r J decreases as r max decreases; see 
figure 0). 
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Figure 12. Predicted X-ray core radii, r^, as a function of the emission- 
weighted mean te mperat ure. Tx, in comparison with the o bserva tional data 
from Mohr' et al. (1999) (circles) and Finoguenov et al. (2001 1 (squares). 
The solid lines represent a = 1, while the dashed lines represent a = 3/2. 
The thin lines are calculated with the predicted X-ray surface brightness 
profiles fitted to the j3 model out to rsnn ■ The thick lines take into account 
the observational selection effect on the maximum radius used for the fitting 
(see figure [n| or equation (p4|)). 



Figure 13. Predicted ratio of overdensity radii, rson , to X-ray core radii, 
fc , as a function of the emission-weighted mean tempe rature , Tx, in com- 
parison with the o bserva tional data from Mom' et al. (199S) (circles) and 
Finoguenov et al. (2001 1 (squares). The solid lines represent a = 1, while 
the dashed lines represent a = 3/2. The thick lines plot rsno /t£ with tak- 
ing into account the observational selection effect on the maximum radius 
used for the fitting (see figure |ll| or equation (p4|)). The thin solid lines show 
the concentration parameter, c = r v i r /r a . 



We find that the predicted core radii from our model agree 
with the observational data. The predicted core radii for a — 3/2 
are systematically smaller than those for a = 1. Qualitatively, this 
is what one would expect, as the former profiles are more centrally 
concentrated than the latter ones at small radii (the rescaling of con- 
centration with 1 .7 only brings the outer profile in agreement, while 
in the centre differences remain). We also see this trend in figure^, 
the predicted X-ray profiles. Both the models (a = 1 and 3/2) give 
acceptable fit to the observations, given the large scatter of the data. 

We have tried other concentration parameters. We find that the 
concentration parameter affects the core radii less than the inner 
slope does; thus, inner gas profiles are a potentially powerful probe 
of dark matter distribution at small distances. Our polytropic as- 
sumption may, however, be no longer valid in the central region, so 
that we need to generalize it to account for the observed decrease of 
the temperature towards the centre (McNamara et al. 2000; Allen et 
al. 2001a, 2001b; Arnaud et al. 2001a, 2001b; Kaastra et al. 2001). 
We plan to investigate this in the future. 

Figure |l2| shows the dependence of the X-ray core radii on 
the emission-weighted mean temperature. The positive correlation 
seen in the figure is due to the hotter haloes being larger. 

To investigate the concentration dependence on the mass or 
temperature, we plot in figure [o] r$oo /rf as a function of the 
emission-weighted mean temperature. Our predictions for a = 1 
agree with the data. There is a marginal tendency for cooler haloes 
to have a higher rson/Vc , and hence a higher concentration. We 
also plot the concentration parameter (thin lines) in the figure; it 
shows the same correlation as r^oo/r^, supporting our interpreta- 
tion. 



We find rsoo/r* ~ 0.4c for a = 1, and rsoo/fc 



0.1c 



for a = 1.5. The former gives r c 

x 



0.2r s , while t he lat ter gives 
r~ c ~ ~ 0.05r s . The former is similar to Makino et al. ( 199S| ) despite 
our model not assuming the isothermality. Note that their predicted 
core radii are systematically smaller than the observational data, 
because the concentration parameter they use is larger than that we 
use on cluster mass scales. 

An important observation from figure [l2| is that the data show 
fairly large scatter in rsoo/r^. While this could be explained to 
some extent by the uncertainties in the observational selection ef- 
fect, it could also imply that the universal gas profile as proposed 
in this paper can only be valid in a statistical sense, and thus indi- 
vidual clusters may deviate from it significant ly. This is consisten t 
with the findings from numerical simulations ( Thomas et al. 2001 ), 



but has not been shown with the observational data on cluster mass 
scales. 



4.4 X-ray outer slope /3x 

Figure [l4| plots the X-ray outer slope, fix, as a function of the 
emission- weighted mean temperature, Tx. The observational data 
show a positive correlation between 0x a nd Tx , which is consistent 
with what Horner, Mushotzky & Scharf (1999) and Lloyd-Davies, 



Ponman & Cannon (200C) find in a differe nt observational cata- 



logue. Similarly, Neumann & Arnaud (199S) find a positive corre- 
lation between /3x and X-ray core radii. Since the X-ray core radius 
increases with Tx (see figure [l^), it follows that fix also increases 
withTx. 

In contrast, our model predicts a negative correlation between 
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Figure 14. Predicted X-ray outer slope, /3x> as a function of the emission- 
weighted mean te mperat ure. Tx, in comparison with the o bserv ational data 
from Mohr et al. (19991 (circles) and Finoguenov et al. (2001). The solid 
lines represent a = 1, while the dashed lines represent a = 3/2. The 
thin lines are calculated with the predicted X-ray surface brightness profile 
fitted to the /3 model out to r'soo- The thick lines take into account the 
observational selection effect on the maximum radius used for the fitting 
(see figure |ll| or equation (p4|)). 



fix and Tx, if we fit the X-ray profile to the f3 profile out to a 
fixed overdensity radius such as rsoo- This is easily understood: 
the smaller the mass is, the more centrally concentrated the profile 
becomes, and the larger /3x we obtain. We must, however, also take 
into account the surface-brightness-cutoff selection effect. Once we 
take it into account, the predicted /3x comes to agree with the data 
very well. We show this explicitly in figure |l4| The selection effect 
makes the maximum radii used for the fit smaller for cooler haloes, 
and thus makes flx smaller (/3x decreases as r max decreases; see 
figure^. 

The selection effect would particularly be prominent for th e 
group sample (T x ~ 0.4 - 2 keV) of Helsdon & Ponman ( |2OO0| ), 
for which r max is much smaller than rspp. They find fix ~ 0.3 for 



the group sample (Helsdon & Ponman 2000), much smaller than 
fix ~ 0.6 for clusters. 

We conclude that the observational selection effect causes the 
positive correlation between fix and Tx- We predict that the corre- 
lation becomes negative, if we fit X-ray profiles to the /3 profile out 
to a fixed overdensity radius. 



SURFACE BRIGHTNESS PROFILES OF THE 
SUNYAEV-ZEL'DOVICH EFFECT 



The Sunyaev-Zel'dovich (SZ) effect (Zel'dovich & Sunyaev 1969) 



is now established as a very powerful observational tool for imag- 
ing clusters (e.g., Carlstrom et al. 2000). Measurement of surface 
brightness profiles of the SZ effect gives projected gas pressure of 
clusters. To extract the gas density profile from the measurement, 
we need to deproject the measured SZ surface brightness profile; 



Figure 15. Predicted surface brightness profiles of the Sunyaev-Zel'dovich 
effect (solid lines), and the best-fitting spherical /3 profiles (dashed and dot- 
ted lines). The fit is performed out to either the virial radius (dashed lines) 
or r5oo (dotted lines). Note that rsoo is half the virial radius (see figure |lj). 
The top panels use the dark matter density profile for a = 3/2, while the 
bottom panels use the one for a = 1. From left to right panels, each panel 
shows the case of M = 10 15 , 10 14 , and 10 13 Mq, respectively. 



thus, we need an appropriate model to parameterize the gas density 
and the gas temperature profiles. 

As often done in X-ray observations, the surface brightness 
profile of the SZ effect has also been parametrized with the (5 pro- 
file, 



Mr) = Jsz(0) 



1 + 



r 

,.sz 



l/2-3/3 sz /2 



(36) 



We compare our predicted SZ profiles with the f3 profile. We cal- 
culate the SZ surface brightness profile from 



hz{x) 



dl P ei 



dl 



J/gas 



+ P 



(37) 



Figure [ij] shows the predicted SZ surface brightness profiles 
together with the best-fitting /3 profiles. The fit is performed out 
to either the virial radius (dashed lines), or rsoo (dotted lines) at 
which the dark matter density is 500 times the critical density of 
the universe. The SZ profiles are substantially shallower than the 
X-ray profiles (compare this figure with figure m. This comes from 
the different weights in the gas density for the SZ and X-ray pro- 
files: the SZ surface brightness is weighted by p gas , while the X-ray 
surface brightness is weighted by /9g as . 

Although the (5 profile gives a reasonable fit to our profile, 
the slope of the fitted profile is quite sensitive to the outermost ra- 
dius used for the fit. This tendency is also seen in fitting the X-ray 
profiles to the j3 profile, indicating that a simple power-law profile 
such as the /3 profile describes neither the X-ray nor the SZ profiles 
properly. 

If the gas is isothermal, then the j3 profile implies that the gas 
density profile is given by 
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Figure 16. Comparison between fitted parameters of the /3 profile obtained 
from the X-ray and the Sunyaev-Zel'dovich surface brightness profiles. The 
ratio of the SZ outer slope to the X-ray outer slope, /3sz//3x. is plotted in 
the top panel, while the ratio of the SZ core radius to the X-ray core radius, 
r c Z / r ?> ' s plotted in the bottom panel. They are shown as a function of 
the maximum radius, r max , out to which the surface brightness profiles are 
fitted to the fi profile. The maximum radius is in units of the virial radius. 
The solid lines, dashed lines, and dotted lines indicate the virial mass of 
10 15 , 10 14 , and 10 13 Mq, respectively. The corresponding concentration 
parameter, c, to the mass is quoted. The dark matter density profile for a = 
1 is used. 



(r) = Pgas(0) 



1 + 



(0 



-3/3/2 



(38) 



with 7v = rt 



c and /3 = /3x = /3sz. As we have shown, our 
predicted profiles deviate from this spherical isothermal j3 model 
appreciably. 

One way to quantify how much the /3 model deviates from our 
model is to calculate Psz/ fix and rf z /r^. If the spherical isother- 
mal (3 model predicted exactly the same surface brightness profiles 
as our model, then these ratios should be exactly 1 . FigureMa shows 
predicted /3sz / Px in the top panel and rf z /r^ in the bottom panel. 
They are plotted as a function of the maximum radius, r max , the 
outermost radius used for the fit. We find that /3sz//3x = 1.15—1.2 
andrf z /r^ = 1.1 — 1.15 for arelevant range of r max ; thus, we pre- 
dict that the SZ profiles give systematically larger /3 and r c t han the 
X-ray profiles (by roughly 20%). Yoshikawa, Itoh & Suto (|l998|) 



also observe this trend for their simulated X-ray and SZ clusters. 
The upcoming SZ observations will test this prediction, and hence 
our universal gas and temperature profiles. 



6 MASS-TEMPERATURE RELATION 

Mass and emission-weighted mean temperature measured for many 
clusters and groups are very tightly correlated. This correlation is 
called the mass-temperature scaling relation. Qualitatively, it arises 
from the virial relation between the virial mass and the virial tem- 



perature, T gas (r v i r ) oc M^/ 3 . Traditionally, workers have used a 
self-similar model with a constant concentration to explain this re- 
lation, which gives the same scaling as for the virial quantities. Hy- 
drodynamic simulations have been used to normalize the relation. 
Our prediction for the mass-temperature scaling relation is 



GfirripM v 

3?*vir 



(39) 



where rjx is the emission-weighted mean normalization factor 
given by equation (^). Figure ^ plots rjx as a function of the virial 
mass. We have fixed the normalization factor at the centre, rj(0), 
with equation (p3|). For comparison with the observational data, 
we use A/500 instead of M V [ T . We calculate A/500 from A/ V i r with 
equation (|l4|). 

Since the concentration parameter in our model depends on 
the mass, the resulting mass-temperature relation is no longer the 
self-similar relation, Tx oc M^. Moreover, the overall normal- 
ization is no longer a free parameter, but fixed, in our model for 
a given universal dark matter density profile. We will show below 
that both features bring our model into better agreement with the 
data than the self-similar model. 

We first begin by discussing the mass-temperature relation for 
hotter clusters than 3 keV. The X-ray observations suggest that 
the observed emission-weighted mean normalization is slightly, 
but systematically, higher than predicted by hydrodynamic simula- 
tions (Horner et al. 1999; Nevalainen, Markevitch, & Forman 2000; 
Finoguenov et al. 2001). It means a higher emission-weighted 
me an tem perature for a given mass. For example, Finoguenov et 
al. (2001 ) find that the sample above 3 keV gives about 50% higher 
nor maliza tion, 7/x(obs) ~ 2.0, than the simulation by Evrard et 
al. (1990 which gives 7/x(sim) ~ 1.3; however, there is scat- 



t er am ong different hydrodynamic simulations: Bryan & No rman 
( |1998[ ) give rjx(sim) = (3/2) /t ~ 1.2, Eke et al. (Jl998[ ) give 

7/ x (sim) = (3/2)/3tm ~ L5 > Yoshik awa e t al. (§0<XJ) give 
7/x(sim) = 7 ~ 1.6, and Thomas et al. (2001) give 77x(sim) ~ 



1.5. Here, the second equality quotes the author's notation. Presum- 
ably this scatter is caused by varying resolution, numerical tech- 
niques, methods of analyzing the simulations, and so on. 

Our model predicts a normalization at the upper end of the 
simulations, rjx ii 1.5—1.8, depending on a and the mass of haloes 
(see figure^)). This means that it fits the observational data reason- 
ably well in the high temperature regime. Figure [n] compares our 
predictions for the mass-temperature relation at z = (solid and 
das hed li nes) with the local observational data from Finoguenov et 
al. (2001). The agreement between the model and the data is rea- 



sonably good, especially for k^Tx > 3 keV, where our model is 
only 20% above the observations. Such small discrepancies could 
easily be explained by either systematic effects in observations or 
by the scatter in the cluster properties, which may increase the av- 
erage polytropic index slightly (we find excellent agreement for 
7 = 1.2). Note that an additional kinetic pressure contribution in 
the hydrostatic eq uilibrium equation, 10 — 20% of the thermal pres- 
sure contribution (Navarro et al. 1995), does not bring our model 



into better agreement; it causes both the data and the model to un- 
derestimate the mass b y the s ame amount. 

Finoguenov et al. ( 2001 ) excised the central region of the clus- 



ters when they measured Tx- Unfortunately, it is not clear what 
inner cut-off radius they have used in the analysis. Hence, we com- 
pare Tx for which the central region inside r cut — 0.01c, 0.05c, 
and 0.1c is excised, with Tx with no excision. 0.1c corresponds to 
340 and 150 kpc for M 50 o = 10 15 and 10 14 Mq, respectively, 
so it is a fairly extreme excision. In reality, the excised region is 
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Figure 17. Predicted mass-temperature scaling relations at z = (solid 
and d ashed lines), compared with the observational data from Finoguenov 
et al. ( 2001 ) (squares). The solid line represents a. = 1, while the dashed 



line represents a = 3/2. 



usuall y sm 
figure |l2|); 



' smaller than the core radius that is about 30 — 300 kpc (see 
; thus, r cu t = 0.1c gives a conservative upper limit on the 
excision effect. We find that r cut = 0.1c reduces Tx by about 5% 
and 10% for M 500 = 10 15 M and 10 14 Mq, respectively. The 
excision affects the smaller haloes more, as they are more centrally 
concentrated. The excision also eliminates the difference between 
a — 1 and 1.5 models. The cut-off radius of below 0.05c has neg- 
ligible effect on Tx compared with the no excision case. Figure |l7| 
plots the predictions with no excision. 

The X-ray observations of cooler (< 3 keV) clusters or groups 
suggest that the slope of the mass-temperature scaling relation de- 

,2/3 



viates fro m the virial scaling relation, Tx oc AIJ^ (N evalainen 
et al. 2000; |Sato et al. 200o( |Finoguenov et al. 2001 )■ It becomes 
steeper for the cooler haloes; the cooler haloes have higher Tx 
than predicted by the virial scaling relation. This is exactly what 
we predict from the model for which the concentration parameter 
decreases with the mass. 

Figure |^ shows how our predicted emission-weighted mass- 
temperature normalization factor depends on the virial mass. The 
trend is in the desired direction: smaller mass haloes have a higher 
normalization, and hence higher temperature than predicted by the 
virial scaling relation. This is due to smaller mass haloes being 
more centrally concentrated than larger mass haloes. 

Figure shows that our predictions for cooler haloes fit the 
observational data better than the predictions of the hydrodynamic 
simulations for which the dependence of the normalization factor 
on the halo mass has not been observed (or studied). The lack of 
resolution in the simulations for the smaller haloes could explain; 
it would not be very surprising, given how sensitive the emission- 
weighted mean temperature is to the central concentration of the 
cluster. We have to verify this with better simulations. 

The predictions with larger a or larger concentration param- 
eters fit the data even better, although in all cases the predictions 



are, on average, still above the data. Within our model, a steeper 
concentration dependence on mass can resolve the discrepancy, 
while an alternative suggestion, the temperature has been raised by 
some energy (or entropy before the infall) injection, is also possi- 
ble (Bialek, Evrard & Mohr 2000). A possible systematic bias in 
the mass determi nation with the hy drostatic equilibrium equation 
could play a role (Evrard et al. 199( ). 



Our results show that the emission-weighted mean tempera- 
ture is very sensitive to the central region of the haloes, and that a 
simple self-similar model cannot be used to compare the observa- 
tions with the models. As we have shown, a more realistic model 
beyond the self-similar model that includes the concentration pa- 
rameter depending on mass removes much of the discrepancy with 
the data. 



7 DISCUSSION AND CONCLUSIONS 

In this paper, we have constructed a model which predicts univer- 
sal gas density and gas temperature profiles from the universal dark 
matter density profile. We assume that the gas density traces the 
dark matter density in the outer parts of haloes. Many hydrody- 
namic simulations have verified the assumption. 

This assumption, together with the hydrostatic equilibrium 
and spherical symmetry, uniquely specifies the gas density and tem- 
perature profiles that are related through a polytropic equation of 
state. The polytropic assumption breaks down in the inner parts 
of clusters (well inside the gas core), as shown by recent Chan- 
dra (McNamara et al. 2000; Allen et al. 2001a, 2001b) and XMM- 
Newton (Arnaud et al. 2001a, 2001b; Kaastra et al. 2001) temper- 
ature profiles. For the outer parts, however, the assumption holds, 
and our predictions based on the assumption agree with both simu- 
lations and observations. 

Our predicted gas density and gas temperature profiles agree 
with the hydrodynamic simulations both in the profile and in the 
overall amplitude. The predicted polytropic index is 7 ~ 1.15, 
which agrees with the observational data from Finoguenov et al. 
(2001) remarkably well. 

We have predicted the core radii and the outer slopes, (3, the 
parameters which empirically characterize surface brightness pro- 
files of X-ray and the Sunyaev-Zel'dovich (SZ) effect. Fitting the 
predicted X-ray profiles to these parameters, we have found that 
our predicted X-ray core radii and X-ray outer slopes, fix, agree 
with the observations, after taking into account the observational 
surface-brightness selection effect. The selection effect makes fix 
increase with the temperature, explaining the observed trend. The 
X-ray core radius is less sensitive to the selection effect, and in- 
creases with the temperature, in agreement with the observations. 

We have calculated the SZ surface brightness profiles, and fit- 
ted them to the /3 profile. The fitted parameters are very sensitive to 
the outermost radius used for the fit, suggesting that the j3 profile 
is not a proper description of the profile. We predict that the core 
radii and /3 fitted from the SZ profiles are systematically larger than 
those fitted from the X-ray profiles, if we use the same r max . More- 
over, for the SZ effect, one may use larger r max than for X-ray in 
the fit; it further increases /?sz compared with fix- 

We have predicted the mass-temperature scaling relation, and 
compared it with the observational data as well as the hydrody- 
namic simulations. We have found good agreement between our 
predictions and the observations, especially for hotter clusters. For 
cooler clusters and groups, we have found that our predictions agree 
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with the observations better than the hydrodynamic simulations, 
which predict much lower temperature than observed. 

The disagreement between our predictions and the simulations 
could be caused by the limited resolution in the simulations, as 
the emission-weighted mean temperature is quite sensitive to the 
gas distribution in the centre. On the other hand, it is also possible 
that some physics is still missing in the simulations, but does not 
show up in our phenomenological description of the temperature- 
density relation. We need careful studies of the numerical resolu- 
tion and ad ditional physics to resolve this issue (S uginohara & Os- 
trike r 1998; pearce et al. 200d; tewis et al. 200c|; Yoshikawa et al. 
2000; [Vluanwong et al. 200 1[ ) 



Some disagreement between our model and the observational 
data remains for the mass dependence of the concentration param- 
eter obtained from pure iV-body simulations. This may be removed 
if a steeper dependence is invoked; for example, the baryonic cool- 
ing contraction and subsequent readjustment of the dark matter 
generically increases concentration, and the effect should be more 
important for smaller mass haloes for which a larger fraction of the 
gas is cooled, and transformed into stars. This trend is in the right 
directi on, and agrees with the observational results of Sato et al. 
fl2000| ) who have measured very steep concentration dependence 
on the mass. We plan to investigate this in more detail in the future. 

We have not attempted to predict other scaling relations 
such as the X-ray luminosity-temperature relation and the central 
entropy-temperature relation. To predict them, we need at least one 
additional uncertain parameter, the gas mass fraction, which is a ra- 
tio of the gas mass to the dark matter mass. Our method does not 
constrain this parameter, as this parameter does not appear in the 
hydrostatic equilibrium equation. While the gas mass fraction may 
be equ al to the cosmic mean baryon fraction, fib/fim, on large 
scales (White et al. 1993), observational evidence suggests that the 



gas mass fraction increases with mass (Arnaud & Evrard 1999 



Mohr et al. 1999 ). This could be explained if some non-negligible 
fraction of the gas were transformed into stars or expelled out of 
the halo. 

Since the X-ray luminosity is proportional to the gas mass 
fraction squared, and the central entropy is proportional to the gas 
mass fraction to —2/3 power, the resultant predictions are signifi- 
cantly affected by the uncertainty in the gas mass fraction. In other 
words, we can explain the luminosity-temperature relation and the 
central entropy-temperature relation without invoking any preheat- 
ing models, if the gas mass fraction varies with the temperature 
(Bryan 2000). Therefore, even without changing our model, we can 
predict these relations consistent with observations, by using an ap- 
propriate gas mass fraction depending on mass. In addition, since 
the size and the density of the core region affect X-ray luminosity 
very strongly, even small modification in the profile leads to large 
change in luminosity; thus, our model does not provide robust pre- 
diction for the X-ray luminosity. 
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